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<N . Abstract 

We present and analyze QCD corrections to the baryon magnetic moments in terms of the one-, 
two-, and three-body operators which appear in the effective field theory developed in our recent 
' papers. The main corrections are extended Thomas- type corrections associated with the confining 

£N| . interactions in the baryon. We investigate the contributions of low-lying angular excitations to the 

o : 

moments quantitatively and show that they are completely negligible. When the QCD corrections 

P-T 

are combined with the non-quark model contributions of the meson loops, we obtain a model 
which describes the moments within a mean deviation of 0.04 /ijv- The nontrivial interplay of the 

i-P ■ 

. two types of corrections to the quark-model moments is analyzed in detail, and explains why the 

• i— i , 

quark model is so successful. In the course of these calculations, we parametrize the general spin 

5-H ' 



structure of the 3 = \ + baryon wave functions in a form which clearly displays the symmetry 
properties and the internal angular momentum content of the wave functions, and allows us to use 
spin-trace methods to calculate the many spin matrix elements which appear in the expressions 
for the moments. This representation may be useful elsewhere. 
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I. INTRODUCTION 



Baryon magnetic moments ( baryon moments for short) have been studied intensively 
for many years using different models and approaches |2], |3|, f|, ||, ||, [7], ||, [K| [11 



12], |T3|, 111, ITR [16], K% 113, [01. It initially appears somewhat surprising that the simple, 



nonrelativistic quark model (QM) [|T], |2|, [| |[] is so good, with more complicated models 
giving only small improvements in the quantitative description of the baryon moments. 
This is now understood. 

In our recent work on the connections between the QM and effective field theory IIR ITU 



we showed that the constituent QM for the baryon moments can be regarded as a rewriting 
of the relativistic chiral effective field theory (EFT) in which only the one-body moment 
operators are retained, a result also obtained by Morpurgo using his general parametrization 
method [H], |17]] derived exactly from QCD. Like Morpurgo, we found that the success 
of the additive QM is due to the numerical dominance of the one-body operators over the 
nonadditive two- and three-body operators, and suggested why the latter operators should, 
in fact, be small. 

In our earlier dynamical model |15j , an initial baryon moment operator was derived from 
quenched QCD using the Wilson-loop approach of Brambilla et al. [p0|| . This appears as 
a sum of single-particle moments /i« ~ (eQi/2Ei), where Ei is the kinetic energy of quark 
i, and a set of Thomas precession terms connected to the binding interactions. In flEfl , we 
studied the further corrections to the baryon moments from from meson loops, which are 
absent in quenched QCD. This combined approach substantially improved the agreement 
between the theory and experiment, leaving an average difference between the theoretical and 
experimental octet moments of only 0.05/ztv. We did not, however, analyze the components 
of the model, or its description in EFT, in detail. The central idea in the context of our 
later work on EFT Jl8|, [19| was to use a dynamical model, namely the QCD-based quark 
model, to obtain reasonable estimates of the input parameters connected to the one- and 
higher- body operators in the effective field theory, then to add the meson loop corrections 
using chiral perturbation theory. 

Our goals here are (i) to present the general structure of the baryon moments in effective 
field theory, parametrized in a form which connects easily to dynamical ideas; (ii) to analyze 
the contributions of meson loops and the dynamical QCD corrections in terms of the EFT; 
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(iii) sketch the relevant parts of our dynamical calculations, including our analysis of the 
contributions of low-lying angular excitations; and (iv) to show how the loop and dynamical 
corrections combine to give a good description of the moments. Items (ii) and (iii) involve 
the structure of the wave functions used to estimate the dynamical corrections, so we give 
their form and the methods of calculation used. Some of this material is potentially useful 
in other contexts. 



The paper is organized as follows: In Sec. [II A| , we briefly review the structure of the 
baryon moments in the effective field theory, where the QM moments arise from one-body 
operators. In Sec. |IIB| we sketch a derivation of the QM moments and further two- and three- 



body dynamical corrections from QCD using a Wilson-loop approach |T5|, jig] . We determine 
the structure of the j = \ + baryon wave functions needed to calculate these corrections 
and develop useful the trace methods for the calculation in Sec. III[ We analyze the 
contributions of meson loops, the dynamical QCD corrections, and internal orbital angular 
momenta numerically in Sec. |I^, and show why both the loop and dynamical corrections are 
needed to obtain a good description of the data. We present our conclusions and discussion 
in Sec. |V|, and give some additional information, such as proofs and sample calculations, in 
the appendices. 



II. THEORY OF THE MOMENTS 

A. Moments in chiral perturbation theory 

In two earlier papers [18|, |19[ , we analyzed the structure of the baryon magnetic moment 
operators to 0(m s ) in heavy-baryon chiral perturbations theory (HBChPT), connecting 
the general spin-flavor structure to the spin dependence of the underlying dynamical theory. 
The analysis was done using flavor-index labeling of the effective baryon fields B^ k (x), where 
i, j, k G u,d, s, and 7 is a Dirac spinor index. The transformation properties of the fields are 
the same as those of the operators 

Bi jk = le abcq r q y k c (c^u (i) 

constructed in terms of free quark operators which carry the color, flavor, and spin structure 
of the baryon. This expression for the relativistic fields reduces in the rest frame of the 
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FIG. 1: One- and two-body contributions to the baryon magnetic moments at the quark level. 
Lines with small wiggles represent photons. There is a factor <x-B <-> —o^F^ in the Lagrangian 
at the quark-photon vertex, where B is the external magnetic field. Crosses correspond to factors 
of the quark charge matrix Q, and dots to insertions of the flavor matrix M. Graphs (a) and (b) 
give the SU(6) structure for the moments. 



heavy baryon to the nonrelativistic structure familiar in the nonrelativistic quark model 

1 



n ijk 



6 



(2) 



and provides a straightforward connection to semirelativistic dynamical models for the 
baryons. It also connects directly to Morpurgo's general parametrization method for spec- 
ifying the most general spin and flavor structure of matrix elements in QCD |7j. In the 
following, we will refer to the structure in terms of "quarks" for convenience, but emphasize 
that the description is completely equivalent to a relativistic effective field theory. 



As shown in ||19H , the octet moments \x can be written through first order in the chiral 
symmetry breaking parameter m s in terms of seven operators m r which correspond to the 




<<\AAs< 



(i) 



FIG. 2: Three-body contributions to the baryon magnetic moments. The two graphs in (i) appear 
together with the same overall coefficient. The structures in graphs (h) and (i) are reducible to 
those in Figs. |l] and |2|(g). 



the diagrams (a)-(g) in Figs. [I] and 0. The contribution of the octet moments to the chiral 
Lagrangian for is given to 0(m s ) as a matrix element of the form (BfiB) ■ B where 



E 



/i r m r 



(3) 
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m„ 



m„ 



is a combination of the quark-level operators 

EWi. ( 4 ) 

i 

£ (QMcr)i, (5) 

i 

£ Qi^ji (6) 

E M <(^)^ ( ? ) 
E ^( M °-)^ ( 8 ) 

m / = E (W^i* (9) 

Ulg = £ W^Jfe- (10) 

Q is the diagonal quark charge matrix, Q = diag(2/3, —1/3, —1/3), and M is the diagonal 
mass matrix for the strange quark, M = diag(0, 0, 1). A flavor index attached to a Q or M 
means that the matrices are taken to act on that flavor quark so that, for example, B(Qcr)iB 
is the matrix element Bk'j'i'Qi'i&j'j&k'k a 'iBijk-i with <Xj acting on quark % in the representation 
in Eq. (fj), — > er^. These operators are shown diagrammatically in Figs. |l]and|2] (Figures 
10 and 11 in ||19|| ). The expressions for the baryon moments in terms of the fi r is given in 



Table ffl \2t 



The coefficients fi r of the operators above are not specified in HBChPT, and the effective 
field theory is not predictive without further input. The additive quark model for the 
moments involves only the one-body operators m a and m^, Figs. [1] (a) and (b). These 
combine to give the effective QM moment operator fiaQ+fibQM. This describes the moments 
quite well, with a root-mean-square deviation of theory from experiment of 0.12 //at, about 
11% of the average magnitudes of the moments. The underlying reasons for this striking 



success have been discussed in |19[] from the point of view of the weak spin dependence 



of the interactions in dynamical models, and in [|7|, [TT], 17| in a QCD-based analysis. These 



related analyses conclude that the two-body contributions to the moments involving m c -mj 
and the three-body contributions proportional to m g — rrij should all be small. Thus, graph 
|l](c) gives a Thomas-type contribution to the moments from the spin dependence of the 
quark interactions as discussed in [1C] and [HJ. Graph |l](d) describes the effect at 0(m s ) of 
the strange-quark mass on the leading contribution to the moment from a different quark, 



TABLE I: The coefficients for the expression of the octet baryon moments in terms of the operators 
corresponding to the diagrams in Figs. H and |, labelled as in the figures, plus the operator niMM = 
Yli^jiQi-MjMk — (QM)iMj] of second order in the strange-quark mass matrix M. The moment of 
the S° is determined by isospin invariance, = ^(Me+ + ME-)- 
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graph |T](a), through the mass dependence of the wave function, while |l](e) and |T](f) give 
the corresponding effects on the Thomas term. Graph (g) in Fig. |2| gives the three-body 
strange-mass correction to the Thomas term. The extra three-body operators (h) and (i) 
arise from the effect of spin-spin interactions on the leading and Thomas terms at 0(m s ). 
These appear in a combined treatment of the octet and decuplet moments, but are reducible 
to the preceding terms if we consider only the octet. 

The dynamical problem in understanding the baryon moments is the actual calculation 
of the coefficients fi r , and of possible terms of 0{m 2 s ) and O(m^) which could affect the 
H, 5*, and Q~ moments. The origin of the leading one-body terms is understood at least 
semi-quantitatively, as discussed below. Two- and three-body contributions to the moments 
can be generated explicitly in HBChPT by meson loop corrections to the leading one-body 
contributions to the moments. The loop corrections have been considered by many authors 
H; i) H Hi 0; HH HI; O; HU • Two- and three-body contributions also appear through the 
nonzero initial values of the coefficients /i c -/U 9 which arise in dynamical models from Thomas- 
type terms associated with short- distance QCD interactions and the long-range confining 



interaction |15|. A model which combines these calculations |16| gives an excellent fit to 
the data, with a mean deviation of the calculated moments from experiment of only 0.05 
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/xjv- In the following sections, we discuss the origin of the dynamical corrections to the QM 
moments, show that possible orbital contributions to the moments are negligible, analyze 
how the loop and dynamical corrections combine to produce the good fit found in |16[], and 
point out where problems remain. 



B. Baryon moments in a QCD-based quark model 



In previous work |T5|, 0, we derived the QM for the baryon moments, including dynam- 
ical corrections, in the context of QCD. Our approach was based on the work of Brambilla 
et al. pOfl , who derived the interaction potential and wave equation for the valence quarks 
in a baryon using a Wilson-line construction in quenched QCD. Their basic idea was to 
construct a Green's function for the propagation of a gauge-invariant combination of quarks 
joined by path ordered Wilson-line factors 

U = Pexp(ig A g ■ dx] , (11) 



where A g is the color gauge field. With internal quark loops omitted (the quenched ap- 
proximation to QCD), the Wilson lines sweep out a three-sheeted world sheet of the form 
shown in Fig. |3] as the quarks move from their initial to their final configurations. The 
approximation ignores the effect of meson loops. 

By making an expansion in powers of 1/rrij using the Foldy-Wouthuysen approximation 
|22|1 , and considering only forward propagation of the quarks in time, Brambilla et al. were 



able to derive a Hamiltonian and Schrodinger equation for the quarks, with an interaction 
which involves an average over the gauge field. That average was performed using the 
minimal surface approximation in which fluctuations in the world sheet are ignored, and the 
geometry is chosen to minimize the total area of the world sheet subject to the motion of the 
quarks. The short- distance QCD interactions were taken into account explicitly. Finally, 
the kinetic terms could be resummed. The result of this construction was an effective 



Hamiltonian [£C| to be used in a semirelativistic Schrodinger equation = E^, 



H = Y. \/p 2 i+ m i + a ( r i + r 2 + r s) - \— + Vsd > ( 12 ) 



- 3 Tij 
Kj J 



where Vsd is the spin-dependent part of the potential. Here = £Cj — Xj is the separation 
of quarks i and j, is the distance of quark i from point at which the sum r\ + r 2 + r 3 
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world sheet 




Wilson line 



FIG. 3: World sheet picture for the structure of a baryon in the Wilson-loop approach. 

is minimized, the parameter a is a "string tension" which specifies the strength of the long 
range confining interaction, and a s is the strong coupling. 

We will make one further approximation, known to be good flST] , and will replace (r\ + 
^2 + ^3) in the confining potential by ^(r 12 + r 23 + r 31 ) and make corresponding changes in 
the associated spin- dependent terms. This approximation allows simple analytic calculation 
of a number of matrix elements we will need. With this change, the Hamiltonian becomes 



H 
H 



SD 



Hq + Vsd, 



i<j 



3r, ; 



a Si • (r 12 x Pl ) - — Sx ■ (r 13 x Pi ; 



I 



Si 



a. 



4mf ri3 
r 12 x Pi)— + (ri 3 x pj — 

'12 '13 

2 1 « s 



3m 2 

2 — ]_ _ x p — . ri3 x p 

3 mirri2 r\ 2 & mim 3 r± 3 



2 a s 



VTLiVTL2 3 r 3 

1 16vr 



mim 2 9 



~2-(Si • r 12 )(<S , 2 • r 12 ) — Si ■ S 2 
1 

c^s^ 3 ( r i2)Si ■ S 2 + permutations + 



12 L r l2 



(13) 
(14) 



(15) 
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with Si = <Tj/2 the spin operator for quark i. The masses which appear in the spin-dependent 
terms are to be interpreted as effective masses, with l/m^ ~ l/Ei, and are not necessarily 
equal to the masses with the same labels that appear in Ei = ypf + mf . The delta function 
spin-spin interactions will not play a role in the following, and will be dropped. The terms 
hidden in the ellipsis are momentum- dependent interactions of 0(p 2 /m 2 ) relative to the 
terms given explicitly. These lead only to small corrections to the already-negligible orbital 
contributions to the moments, and will be ignored. This model has been used in detailed 
calculations in and [Q to obtain good fits to the baryon spectrum up to ~ 3 GeV. 

To derive a expression for the baryon moments in the same approximation to quenched 
QCD, we redo the calculation of Brambilla et al. with the gauge interaction extended to 
include the electromagnetic vector potential A em (xi) = B x x q /2 associated with a constant 
external magnetic field B. We can then pick out the modified magnetic moment operator 
through the relation 

AH = -fiB. (16) 

The result is 

M=E(^ QM) +^r)+^- (it) 

j 

Here fij is the QM moment operator fij&j, while A/j^ M and fii are the leading cor- 
rections to the baryon moments associated with the binding interactions in VsD and the 
nonzero orbital angular momentum in the baryon, respectively. The latter two can be ob- 
tained directly by making the minimal substitution p { — ► — ejA(a3j) in Eq. (|T2| ) and 
isolating the B-dependent terms. This gives 



Afi QM = 111 *. 

6mi 


[(«!■ 


• r 12 ; 


)<Tl- 


(xi ■ 


■ <Tl, 


)r 12 ] 


+ (2 


«-> 3) 




2mi rf 2 




• r 21 ; 


)<Ti- 


(x 2 


■ <Tl, 


)r 2 i] 


+ (2 


<-> 3) 




A*i o" 
4mi ri2 


[(*!■ 


■r 12 ] 


)(Ti- 


(xi 


■<*!. 


)ria] 


+ (2 


^3), 


(18) 



where \ii = tijlmi with m« an effective mass. The results for A/j,^ and A/X3 which can 
be obtained by cyclic permutation of the indices. Finally, 

t*L = ^2lijX j xp j = L j (19) 

j 

where Lj is the orbital angular momentum of the j-th quark. 
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As discussed in | 19| , A/j,® m arises from a set of Thomas precession terms. The diag- 



onal terms in the first and third lines of Eq. ( |T8| ) are associated with the spin-same-orbit 
interaction in Eq. (|15|), and give a largely multiplicative correction to /ij. The two-body spin- 
other-orbit terms introduce new, nonadditive structure, and are important in improving the 
simple quark- model fits to the moments [|16|| . These terms are proportional to the short- 
distance spin-dependent part of the potential of order a s divided by EiEj, so are expected 
to be small. 

To obtain an approximate operator form for the QM moment fij to use in later cal- 
culations, we act on the electromagnetic interaction term —eon ■ A{x,\) with the Foldy- 
Wouthuysen transformation which reduces the operator Dirac Hamiltonian (a • p + l3m) for 



the free motion of a quark to the two-component form H = f3E [p2 |, with E = \Jp 2 + m 2 the 



kinetic energy operator in Eq. (fl^) . This leads to 0(e) to a two-component spin-dependent 
interaction 

:cr B 6 / cr ■ (B x p) x p. (20) 



2yJPj + m | 4^p 2 + m 2 (^p 2 + m 2 + rrijj 

The first term can be identified with the QM moment interaction —fijcr ■ B and gives 

N = (e j /(2E j )) (21) 

when averaged over the momentum distribution in a baryon. We will adopt this form in 
later calculations to obtain reasonable estimates of the effect of the different binding of the 
quarks in different baryons on the fij, but will otherwise treat the fij as parameters. In this 
treatment, the second term in Eq. (|20|), which behaves similarly as far as the averages are 
concerned, gives a single-particle contribution to the baryon moments which is absorbed in 
the adjustment of the parameters fij ~ ej/2rrij in terms of the effective masses. 



III. BARYON WAVE FUNCTIONS 



A. Spin structure of the baryon wave functions 

To estimate the dynamical corrections to the baryon moments, we need approximate wave 
functions for the Hamiltonian in Eq. ([12]) . We have found it very useful to construct the wave 
functions for the j = \ + baryons in a form that allows us to use trace methods to calculate 
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the many spin matrix elements which appear in the expressions for the magnetic moments. 
This construction also has the advantage of displaying clearly the symmetry properties and 
the internal angular momentum content of the wave functions. The basic idea is simple. We 
start with the j — \ spin wave functions for three quarks in the quark-model ground state, 
and construct all further wave functions by the application of even-parity scalar operators 
with the correct symmetry to give allowed states. Because the total angular momentum 
operator J commutes with scalar operators, the new wave functions will retain the original 
eigenvalues j(j + 1), j 3 of J 2 and J 3 . Our starting spin wave functions are of the familiar 
form 

x s A =1 = ^=[(T! + IT)T-2(TT)!], 

2>2 vb 

x? 12 * 1 = -^=[(U + IT)I-2(U)T], 

2' 2 

xf 12 r° = 4=(TWT)T, 

2'2 V 2 

x f 12 =° = J^ti-ITH, (22) 

2' 2 y 2 

where we will label the quarks 1, 2, 3 in order. The spin functions for S\2 = 0, 1 are 
connected by the operator (<Xi — <x 2 ) ■ <x 3 which changes the symmetry in spins 1 and 2, 

as is easily checked by direct calculation. 

The functions xf 1 ^ ' 1 & ve a complete basis for the j = \ s P m states for three quarks. A 
general j — \ wave function can therefore be written in the form 



xl^^^xl 1 ]; 1 (24) 



where <pi and 02 are scalar operators symmetric in <j\ and <r 2 . We will construct the full 
operator directly; the part antisymmetric in ct\ and cr 2 will generate the = component 
of ip when acting on the symmetrical Su = 1 spin function. 

We will concentrate initially on baryons containing two like quarks, taken as 1 and 2 in our 
labeling. The singlet color wave function of the baryon is antisymmetric in the interchange 
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of any two labels, so the space-spin wave function and, therefore, the operator <\> must be 
symmetric under the interchange of 1 and 2. No symmetry requirement exists with respect 
to the third quark. We will work in the center-of-mass system of the three quarks, and use 
the Jacobi coordinates p and A, defined respectively as the separation of quarks 1 and 2, 
and the separation of quark 3 from the center of mass of 1 and 2, to describe the internal 
configuration of the baryon. It is also useful to introduce Jacobi coordinates in which the 
pairs 2,3 and 3,1 are singled out, with 

P = r i2) A = r 12)3) 

p' = r 23 , A' = r 23 ,i, (25) 

P" — r 31; A" = r 31;2 , 

where = Zj — Zj, = Ry — Zk, and = (rriiZi + mjZj)/mij. The relations among 
these coordinates and among the corresponding momenta are given in Appendix [AJ We note 
here that p and A are, respectively, odd and even under the the interchange of 1 and 2. 

It will be useful before determining to determine the possible orbital angular momentum 
configurations for a spin j = | + baryon. We introduce two angular momenta, L p between 
quarks 1 and 2, and L\ between quark 3 and the center of mass of 1 and 2, and combine 
these to produce the total orbital angular momentum L. This must combine in turn with 
the total quark spin S = |, | to give j = | + . Clearly, L < 2, while L p and L\ must both 
be even or both odd to give a positive parity state. If we restrict our attention to states 
with L p , L\, L p + L\ < 2 as well, the complete wave function can be constructed using the 
configurations 

(L„M = ( W '» )(M ' » (1 ' M) (26) 
{ (1,1,2) (2,0,2) (0,2,2) 

Each unit of L p (L\) in the wave function corresponds to one factor of p (A) in the corre- 
sponding angular momentum tensor, and a symmetry factor -1 (+1) for an interchange of 
quarks 1 and 2. The required symmetry of the wave function under quark exchange then 
limits 5*12 to the values S12 = 1 (S12 = 0) for L p even (odd) for the baryons with quarks 1 
and 2 identical. This eliminates the configuration (1,1,2), S12 = since this combination of 
angular momenta gives j = 3/2, 5/2 only. 

It is now straightforward to construct the operator in Eq. fl24|) using the Pauli spin 
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matrices <Xj, the coordinate vectors p, A, and the tensor tij, 

Uj(x) = 3cTj X (Tj X — (Ti-CTj x 2 . (27) 

The possible components are given in the labelling (L p , L\, L)S\ 2 by 

(0,0,0)1: 1 (28) 

(1.1.0) 0: (<n - <r 2 ) ■ <7 3 p ■ A (29) 

(1.1.1) 0: (ffi-ff 2 ).pxA, (<Ti - a 2 ) x <r 3 • (p x A) (30) 
(2,0,2)1: t 12 (p), t 23 {p)+t 31 {p) (31) 
(0,2,2)1: t 12 (A), t 23 (A)+t 3 i(A) (32) 



Two operators are shown in the rows for L = 1, 2. However, it is clear that the two cannot 
be independent for the j = | + baryons because there is only one way to reach j ; = | by 
combining the indicated values of L, Si 2 , and s 3 = \. We show directly in Appendix [FJ 
that the second operators are not independent of the first for j = | (this is not the case for 
L = 2, j = |). We will therefore write in terms of the first operator in each row, with 

= a + ib (<Ji - <T 2 ) ■ p x A + ct 12 (p) + dti 2 (X) + e(cr 1 - <r 2 ) ■ er 3 p ■ A, (33) 

where a, . . . , e are functions of p 2 and A 2 only P7fl . In the general case, the coefficient 
functions can also depend on (p • A) 2 , but this requires the introduction of values of L p = 
L x >2. 

The analysis above holds for quarks 1 and 2 identical, that is, for the p, n, S 1 * 1 , 5°, and 
S~ baryons. The S° and A are both uds systems, and are distinguished by the interchange 
symmetry of the ud pair. This is even for the E°, so the S° wave function has the structure 
given in Eq. (0). The A is different, with odd ud interchange symmetry. An analysis similar 
to that above gives the spin dependence of the A wave function as ipf . = 0Axf 12- \ with 

0A = a A (T 3 ■ (<T! - <T 2 ) + b A p ■ A + C A i(T 3 ■ p X A 

+d A p x A • <r 3 x (o-i + <r 2 ) + e A t l2 {p, A), (34) 

where 

ti 2 (p, A) = i[3cr r pcr 2 -A- cr 1 -cr 2 p-A+ (1 <-> 2)]. (35) 
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B. Trace methods 



Since the total angular momentum operator J commutes with the Hamiltonian H, the 
matrix elements of H satisfy the relations 



(j'J' 3 \H\jJ 3 ) = U'-.h 11 J-J-^f/j^r ( 36 ) 
(j,j' 3 \H\j,j' 3 ) = (j,h\H\j,j 3 ), (37) 



so, for j = |, j 3 = ±\ 



^k\ H %h) = -[(\,\\H\\,\) + {\,-\\H\\,-\)\ 

= h^i-lB- (38) 

A J — 2 

We can use the projection operators to get rid of the restriction to j — \ in the trace. 
Formally, 

Ps,s 12 = ^2\Si2;sJ 3 )(S 12 ;s,j 3 \ . (39) 

33 

Now P^iXj 3 2=1 = Xj.j 2=1 , while Pi x annihilates ^|j 3 , so by introducing the projection 
operator 

P| ,i = ^ [3 + <T! • <r 2 - 2<n ■ <r 3 - 2<r 2 • <r 3 ] (40) 

in matrix elements, we can use a basis of independent states j, j for quarks 1,2, and 3 and 
sum over the spin indices without restriction. This gives 



= 5 E ftWPfr 



2 

Al,A2,A3=±2 

= l -Tv^H4>P hl (41) 
where Tr = Tr!Tr 2 Tr 3 and \ is the spin projection of the i-th quark (i = 1, 2, 3). Note that 

Trl = 8, 

^TrPi ;1 = 1, (42) 

Tr (Ti = Trcr 2 = Tr<r 3 = 0. 
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C. Approximate wavefunctions 



To get approximate wave functions to use in our calculations of the binding corrections to 
the moments, we have made variational calculations of the energies and approximate wave 
functions of the ground state baryons and their first excited states using the Hamiltonian 
in Eqs. (|T2|) and (|1^). Our calculation and results are summarized below. Some details are 
given in Appendix A. 

We emphasize that we are not trying to make detailed, accurate calculations, but rather 
to obtain crude wave functions which allow us to estimate the small binding corrections. 
More accurate calculations would probably only be useful in conjunction with further im- 
provements in the underlying theory. 

We rewrite the baryon wave function in Eq. fl33D as 



^§, m = W + bo^Pb + c ^ c + dafa + eoi)Xi^ , (43) 
where ip^ is a normalized wave function of the ground state and 

?/> b = lN b (<7i - (T 2 ) ■ p X XljjQ 

4 = N c t 12 (p)r (44) 

h = N d t 12 (X)^ 

= N e ((Ti - cr 2 ) ■ cr 3 p ■ A V>o , 

are the normalized wave functions of the excited states. Now if we choose ip^ , ■ • • , ip^ as the 



Gaussians 28 



Vtf = ( ^ ) exp [-{flf p 2 + pfx")} , A = a, ■ ■ • , e , (45) 



3/2 



then the normalization coeffcients Nf, , ■ ■ ■ ,N e are 

Vf' Nc ~W^ d ~73T Ne ~ — - (46) 

The coefficients bo , • • • , eo can be evaluated using perturbation theory in which the excited 
states are considered as the perturbation states. In such a way, we have 

= A = i,-,e, (47) 
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where V$d is the spin- dependent potential in Eq. ([15|) with the spin-spin and higher order 
terms omitted, and 

E A = $ A \H Q $ A ) , Eo = mH \M). (48) 



H is given in Eq. ([U 

To evaluate the coefficients bo ,■ ■ ■ , eo numerically, we first minimize the energies E A and 
Eq by varying the parameters (3 P , j3\, and thus obtain the best variational wave functions 
ipQ, tjj e of Gaussian form. The calculation uses a s = 0.39, a = 0.18GeV 2 , m u = 

rrid = 0.200 GeV, and m s = 0.500 GeV. The resulting values of (3 P and (3\ are given in Table 
These values differ for different baryons, a point which will be important later. The 
values of excitation energies Ea — E obtained this way and shown in Table [TT] are generally 



consistent with those found in the more elaborate calculations [^] [f29[ . 

These results together with those from our calculation of the matrix elements (^aIKsdIV'o) 
between the ground state and lowest excited states b , • • • , e allow us to estimate the coeffi- 



cients b , • • • , e . We list the values found in Table [TV. 



TABLE II: The energy differences AEa in MeV between the baryon ground states ipo and the 
excited states ipA, A = b, . . . , e in Eq. (f45|). 



Baryon AE b AE C AE d AE e 

N 870 821 820 789 

S 841 837 750 755 

E 810 714 806 738 

n 779 732 730 701 



While the Hamiltonian does mix orbitally excited states into the L p = L\ = L = 
quark-model ground state, the coefficients are very small, ranging from essentially zero to 
about 0.02 depending on the baryon. Because these coefficients only appear quadratically in 
the baryon moments, the orbital contributions to the moments will be completely negligible. 
This result is consistent with the Lichtenberg's finding for the nonrelativistic QM p4 |. 
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TABLE III: Values of the parameters (3 p and (3\ in the Gaussian approximation to the ground state 
wave functions of the octet baryons, Eq. (|^), calculated for m u = ma = 0.200 GeV, m s = 0.500 
GeV, a s = 0.39, and a = 0.18 GeV 2 . 



Baryon 


N 


s 




n 


(3 P 


0.340 


0.347 


0.387 


0.394 


/3a 


0.393 


0.424 


0.420 


0.455 


TABLE IV: The 


coefficients bo, ■ ■ ■ , eo are 


evaluated using 


perturbation theory for a s 


= 0.39, 


a = 0.18GeV 2 , m u 


= m d = 0.343 GeV, and 


m s = 0.539 GeV. 






Baryon 


bo 


CO 


d 


eo 


N 





-0.013 


0.013 


-0.021 


E 


0.012 


-0.013 


0.011 


-0.016 




-0.006 


-0.009 


0.011 


-0.020 


n 





-0.009 


0.010 


-0.016 



IV. ANALYSIS OF THE BARYON MOMENTS 
A. The leading approximation to the moments 

The leading approximation to the QM moments in Eq. (|2~l|) gives fij = (ej/Ej) = 
(ej/ + m'j)). The matrix elements (1/Ej), customarily written in terms of an effective 
mass as ej/rrij, will actually differ for a given quark depending on the baryon in which it 
is confined. To estimate this effect, we have calculated the relevant matrix elements using 
the Gaussian variational wave functions and the methods sketched in Appendix 0. The 
calculation was done using m u = = 0.200 GeV and m s = 0.500 GeV in H , values which 
give both a good ground-state baryon spectrum and reasonable values for the moments. We 
then estimate the corrections to fi u by taking the ratios of the calculated matrix elements in 
the E + , A, and S° to the matrix element in the proton, and multiplying by the actual value 
fi u = 2fi a /3 determined from the data. The correction to fi s in the E states is determined 
similarly using the £ as the reference state. The corrections are small, with, for example, 
[l u decreasing by 0.046 in the S and 0.088 fi^ in the 5 relative to its value in N. 
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The overall corrections to the parameters fi a , ■ ■ • , Hmm are given in Table |Vl[ The largest 
correction is to fid- This arises, as suggested by Fig. |l](d) from the change in the moment of 
one quark, relative to the values for the rrii given above, which arises from a change in the 
wave function when a second quark is strange. 

The nonzero value of [lmm found in this calculation corresponds to the introduction of 
the new structure 

™-B,MM = [QiMjM k - li(QM)jMk} a- (49) 

at the baryon level, where a is twice the baryon spin operator. This is a mixed two- and 
three-body operator. The operator coefficient of cr vanishes for all states which do not 
contain two strange quarks, and has the value 1 for the S° and for the 5~. There are no 
further independent structures for the octet baryons . 



B. Binding corrections 

Without the corrections from the internal orbital angular moments, the moment of a 
baryon B is given by 

»B = + A ^f)(^)B = l4 M + J2 A »!(^*)b, (50) 

3 3 

where the sum is over the quarks in the baryon and we have quantized along B, taken along 
the z axis. The spin expectation values are to be calculated in the baryon ground state. The 
correction AfiJ to the moment of quark j depends on the baryon B in which it appears. The 
final baryon moments depart from the quark model pattern only when the ratios Afi? 
differ in different baryons. 

The general result for the moment operator given above can be simplified considerably 
for the L = ground state baryons. The absence of any Pauli matrices in the wave functions 
allows us to reduce the operators in A/if M in Eq. (|i"8|) to the components proportional to 
<Tj in calculating the ground-state matrix elements, 

(Arf M )s = (A»? M a z ) B , Afif M = l -Tr a t ■ , (51) 

and the contribution to the quark moment //j in baryon B is just A/if = {Anf M ) B . We 
will write A/if as the sum of terms arising from diagonal and Thomas-type corrections, 

A/if = A/if + A/if. (52) 
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These are given by 



Alibi ■ / ■ e i • / ■ i'H 

where the e's and S's are the ground state matrix elements 

2o^,s / Fi-i ' Xj \ (J I Y^j ' X{ 

Two quarks always have the same mass in each octet or decuplet baryon, so appear 
symmetrically in the spatial part of a flavor-independent wave function. We will label these 
quarks 1 and 2, let m\ = rri2 = m, and take 3 as the odd-mass quark if there is one. With 
this labeling, the spatial matrix elements can be reduced to the small set 

e = ei2 = 621, e' = ei3 = ^23, e = e 3i — ^32, (55) 
£ = S 12 = E 21) £ = S13 = S23, S = £31 = £32, (56) 

and the diagonal correction terms become 

A "? = £<^'- E - E '>' 

A^=^-(s-g). (57) 



m 3 

The Thomas-type corrections are somewhat more complicated 



ei \m 777.3 

A/J2 = »(*e+2Li 
C2 \m 7773 

A/73 = **±^. (58) 
e 3 777 v ' 

Note that, for all masses equal, e = e' = e, and £ = £' = £. The matrix elements differ 
only because of effects of the odd-quark mass on the wave functions. 

We have evaluated the radial matrix elements above using the Gaussian wave functions 
obtained in our variational calculation of the ground state energies for the Hamiltonian 
in Eq. ( |T2j) and the methods sketched in Appendix 0. The results are given in Table [V| 
for a s = 0.39 and a = 0.18 GeV 2 , values taken from fits to the baryon spectrum [^l|, |23| 
using the same Hamiltonian. As is clear from Table [V| the e's and the S's are all similar in 
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TABLE V: The values in GeV of the matrix elements e and S defined in Eqs. ( j55| ) and ( |56| ) 
The matrix elements were evaluated for a s = 0.39, a = 0.183 GeV 2 , m u = rrid = 0.200 GeV, and 
m s = 0.500 GeV for the quark moments fi u = 1.943, fi s = —0.640 obtained from a fit to the data. 



Baryon 


€ 


e' 


e 


S 


£' 


£ 


N 


0.017 


0.017 


0.017 


0.050 


0.050 


0.050 


£ 


0.017 


0.024 


0.011 


0.049 


0.064 


0.031 




0.019 


0.013 


0.023 


0.044 


0.033 


0.060 


n 


0.019 


0.019 


0.019 


0.044 


0.044 


0.044 



magnitude even for nil ^ 1113, so the main effect of the diagonal corrections Afi^ in Eq. 
is a uniform shift in the input values of the fi^s. This is absorbed when fi u = — 2 fid and 
fi s are used as parameters in fitting the observed moments. The only nontrivial two- or 
three-body effect of these terms arises from the differences between the matrix elements in 
the different baryons. In contrast, the Thomas- type corrections AfiJ in Eq. ( |5^| ) are mainly 
two-body contributions and remain nontrivial even when the e's are similar because of the 
different charges of the different quarks. 

In Table [Vl| we shown the values of the parameters fi a - fi g , fiMM f° r the diagonal and 
Thomas-type corrections obtained using the input moments fi u = 1.94 fi N and fi s = —0.64 
fix found in a fit to the data which includes these and the other corrections. The effective 



masses used in Eqs. ([)?]) and (|58|) were defined in terms of the input moments since the 
original coefficients in Eq. (|1^) are of the form Ci/mf or ei/rriimj. The smallness of the 
parameters fi c - fi g for the diagonal corrections compared to the values of fi a and fib shows 
that the Afif are indeed mostly contributions from the one- body operators. In contrast, 
the AfiJ contribute mostly to the two-body coefficients fi c - fi g with the largest contribution 
for the Thomas terms fi c and fif, Fig. |I|. The value of fi g , which is associated with the three- 
body operators, is very small is both cases. Finally, the nonzero values of fiuu obtained 
in both cases arise indirectly from the dependence of the matrix elements in Eqs. (|55D and 



(fidj) on the strange-quark content of the state. 
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TABLE VI: The parameters /i a - fi g and Hmm f° r the estimated binding corrections to the moments, 
the diagonal and Thomas-type corrections defined in Eqs. ( |57|) and (|58|), and for the meson loop 
corrections. The last lines give the total correction and the values of the parameters obtained in a 
least squares fit to the measured moments, and the differences with the experimental uncertainties. 



Type 


Ma 


Uh 
to 


lie 


t u 


/<(■ 




a,, 

r"9 


1^ MM 


A/ 





0.011 





-0.067 


0.011 


0.003 


0.003 


-0.019 


A/i D 


-0.306 


0.231 


0.000 


-0.023 


0.001 


0.000 


-0.008 


-0.012 


A/i T 





0.000 


0.151 


0.000 


-0.006 


-0.082 


0.004 


-0.013 


meson loops 


0.123 


-0.208 


-0.339 


-0.022 


0.143 


0.471 


0.185 


0.010 


A^total 


-0.183 


0.034 


-0.188 


-0.112 


0.149 


0.392 


0.184 


-0.034 


Fit to data 


2.793 


-0.933 


-0.077 


-0.037 


0.098 


0.438 


0.044 


221 






±0.012 




±0.070 


±0.010 


±0.036 


±0.070 


±0.140 


Fit - AyUtotal 


2.976 


-0.899 


0.111 


-0.075 


-0.051 


0.046 


-0.100 


0.255 






±0.012 




±0.070 


±0.010 


±0.036 


±0.070 


±0.140 



C. Contributions of internal orbital angular momenta 

As already remarked, the contributions of orbital angular momenta to the moments are 
negligible for the ground-state octet baryons. We will illustrate this in the special case in 
which mi = irt2 and fi\ = ^ = m - is given in this case by 

ii L = \i x L p + — (m 3 /ii + 2m/i 3 ) L x , (59) 

where L p = p x p p , L\ = X x p A , and M = 1m\ + m^. 

Using the wave function given in Eq. ( [4*3D and choosing the Gaussian spatial wave func- 
tions discussed above, we calculate (/x) and obtain for the moment of a baryon B of this 
type 

»b = - (fa + A/if ) - - (/i 3 + A/if) . (60) 
A/if and A/if are the contributions from L ^ configurations in the baryon wave function 
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(Eq. ©), 

M^ ( M + M) , (61) 

with 7?. an overlap integral of order unity. The results are quadratic in the small coefficients 
b , . . . , e in Table [IV| and are negligible on the scale of the other contributions. 



D. Meson loop corrections 

As noted earlier, the meson loop corrections to the baryon moments have been considered 
by many authors |, |, |, |, |, § 0, |, g |10|, [TT|, 0, [13|, 0, |15|, 0, ||, || in the context of 



chiral perturbation theory. We will use the loop corrections obtained in JT6[ . In that paper, 
the one-loop corrections were calculated using a meson-nucleon form factor to control the 
high-energy behavior of the loop integrals, which is not reliable in the context of a strict 
expansion in the chiral parameter m s . The results we obtained after correcting an error in 
P~6|] are given in Table [V| as decomposed in terms of the parameters /i a , . . . , /z 9 defined in 



Eq. (^). Their connection with the baryon moments is given in Table [p. More detail on the 
calculations is given in JTBJ . 

The meson-loop corrections are distributed over all the parameters \x a — f^ g , corresponding 
to contributions from the one-, two-, and three- body operators. The largest parameters, 
fi c = —0.339 and /if = 0.471, arise from the two-body Thomas-type contributions in Figs. 
[T] c and |l]f , the second a strange- mass correction to the first. There is also a significant 
three-body contribution to the moments through Fig. |g. 

E. Analysis 

The corrections to the moment parameters for the baryons discussed here are summa- 



rized in Table 0. The final rows in the table give the total calculated corrections to 
Ha, . . . , fig, jiMMi the values of those parameters obtained by fitting the eight measured 
moments exactly, and their difference. The correspondence of the the individual baryon mo- 
ments and the E°A transition moment to the parameters is given in Table |. The differences 
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TABLE VII: The fit to the baryon magnetic moments obtained using the calculated dynamical 



contributions to the chiral parameters /i a , . . . , fx mm given in Table |V|, with fi u and fi s as free 
parameters. 



Baryon 


V 


n 


£+ 


S" 


5° 




A 


S°A 


Calculated 


2.744 


-1.955 


2.461 


-1.069 


-1.278 


-0.598 


-0.607 


±1.522 


Data 


2.793 


-1.913 


2.458 


-1.160 


-1.250 


-0.651 


-0.613 


±1.610 








±0.010 


±0.025 


±0.014 


±0.003 


±0.004 


±0.080 



for fi a and the combination fi a + fi b give the values of 3fi u /2 and — 3/i s , the input quark 
moments given fid = —fi u /2. 

As was noted in fits to the octet moments using fi u and fi s , or equivalently, fi a and 
fib as adjustable parameters give only slight improvements relative to the additive quark 
model when only the binding corrections, or only the meson loop corrections are included. 
Either set of corrections alone has an incorrect pattern of signs or magnitudes relative to 
the fitted parameters. However, when combined, there are significant cancellations, and the 



total result is much closer to experiment. This may be seen from Table [VII| where we give 
the best-fit results for the seven accurately known moments obtained using the calculated 
values of A/z tota i, and allowing fi a and fib to vary. The fit is quite good, with a mean deviation 
from experiment of 0.04 fi^ for fi p , . . . , fi^. The poorly known E°A transition moment has 
very little weight in a complete fit, and is left as a prediction, /i£°A = 1.522 fi^ compared 
to the measured value 1.61 ± 0.08 fi^- 



The advantage of the decomposition in Table |VI] is the possible hints it provides for 
improving the theory. In particular, the most striking deviation of theory from experiment 
is in the value of the leading Thomas-type term /i c corresponding to Fig. |l|(c). This is well 
determined experimentally, with fi n = — |(/i a — aO- The QM ratio fi n j fi v ~ —2/3 is attained 
only for fi c small, a result which is "accidental" in a general parametrization of the moments 
|T7| , p5fl . Here fi c is small because of a cancellation between the meson loop corrections at 
one loop and the Thomas-type correction Afi T . Some further input will clearly be needed 
to obtain better agreement between theory and experiment. 

The calculated value of the related coefficient fi e , is also significantly different from ex- 
periment. The uncertainties in fi d , fi g , and Hmm are dominated by the relatively large 
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uncertainty in the S°A transition moment. The results nevertheless suggest that the three- 
body contribution fi g may be overestimated, and the calculated two-strange quark term 
A 4 mm is too small and of the wrong sign. We note that baryon mass insertions in the meson 
loop diagrams, so far not calculated, will affect these contributions ]TT)|. 



V. SUMMARY 

We have presented and analyzed the QCD corrections to the baryon magnetic moments 
in terms of the effective field theory description in terms of "quark" operators m a , . . . , iximm 
to be used in baryon matrix elements of the form \xb = {BfiB), with Bj- k the effective field 
operator for the baryon in Eq. (|l|). We have then used the connection of this quark picture to 
dynamical models |7], [11], [18], |19|] to estimate the input values of the unknown parameters in 
the chiral expansion using a semirelativistic dynamical model. This model can be derived in 
a quenched approximation to QCD [^DJ and gives a good description of the baryon spectrum 
up to ~ 3 GeV [[H], 23]. This allows us to estimate the changes in the leading QM moments 



in different baryons and binding corrections to the moments identified earlier [TjJ, and to 
study the contributions of orbital angular excitations. 

We find that the binding corrections to the simple QM picture of the moments involve, as 
expected, mainly contributions from the one- and two- body operators, and are important. 
The changes in the QM moments in different baryons are significant but small. The orbital 
contributions to the moments are completely negligible since the terms in the wave functions 
which involve higher angular momenta are very small. The chiral meson loop contributions 
are large, but tend to cancel with the corresponding parameters for the binding corrections 
leading to a fairly successful overall fit to the moment data. Fits using the loop corrections 
alone are much less satisfactory. The deviations of the calculated parameters m a , . . . , m MM 
from those found in an exact fit to the baryons moments suggests where further input will 
be important. 

We note finally that the procedure used in our construction of the j = \ + baryon wave 
functions displays their symmetry properties and the internal angular momentum content 
very simply, and may be useful in other contexts as may the trace methods used to calculate 
the many spin matrix elements which appear in the expressions for the magnetic moments. 
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APPENDIX A: JACOBI COORDINATES AND RELATIONS 

With i 7^ j 7^ k and i, j, and k run from 1 to 3, we define the Jacobi-type coordinate 
systems as follows: 



1. The space coordinates 



Hj — 

m.A 



m,i(xi -x k ) + rrij(Xj - x k ) 



rriij 



rnijUij + m k x k 



R ijk = m ■ (Al) 

Here rriij = m i + ™>j, M = rrii + rrij + m k , and Hij k is the usual center-of-mass coordinate. 
With these definitions, the roles of i, j, and k are completely symmetric. 

Thus, if one defines coordinates p = r 12 , A = r 12j 3, p' = r 2 3, A' = r 2 3 i i, p" = r 31 , and 
A" = r 3 i ;2 then p', A' and p", \" can be expressed in terms of p and A and conversely. For 
example, if m\ = m 2 , we find that 

P =A-f, \' = -—(m 3 \+^- p] , (A2) 
2 m 23 V 2 J 

p" = -A - f , A" = - fm 3 A - . (A3) 

2 m 3 i V 2 J 

One can therefore work with any of the pairs p, A or p', X' or p", A", and switch between 

them as necessary. The spatial volume element is simply d 3 Rd 3 pd 3 X, or equivalently for the 

other pairs of internal coordinates. 
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2. The momentum coordinates 



Pij ~ i "ij ~ Pi i Pj > 

mfcPjj - m ijPk 

Pij,k - , r ijk - fij + p k , (A4J 

where Py*. = P is the total momentum. 

Now, if one denotes p p = p 12 , p A = Pi 2 ,3> P P > = P23, Pv = Pa3,i> P P " = Psi> and 
Pa" = P312 then one can choose to work with either p p , p A or p p x , or p p „, p x „ since 
there exist relations between them. For example, when mi = m 2 = m the coordinates p p /, 
p A / and p „, p A „ can be expresses in terms of p p , p A as 

Pp' Pa Pp „ _ „ Pa 

^ ~ m~ x 2m~/ Pa '~ Pp T' 

JV =_Pa__^, PA-Pp-f , (A5) 
rrip" m\ 2m p F 2 

where m p = m 1 m 2 /m 12 = 2/m, m A = m 3 m 1 2/m 12 3, m p i = m 2 m 3 /m 2 3, and m p » = 
m 3 mi/m 3 i. The volume element in momentum space is d 3 P d 3 p p d 3 p x , and equivalently 
for the other pairs of internal momenta. 



APPENDIX B: TRACE RELATIONS AND INDEPENDENT OPERATORS 

We present here some trace relations derived from the trace method given in subsection 



[JIB, These results are useful to an evaluation of the spin matrix elements. We also use 



them to prove the equivalence of some states shown in Eq.( 
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1. Results for spin traces 



Let A, B, and C be any spin-independent vectors and x, y, and z be the coordinate-type 
vectors (such as p, p', A, ...), we have 







Tr( CTl - 


- c 2 ) ■ cr 3 crx ■ 


""2 Pi ! 

2 1 


= 0, 


(Bl) 






Tr (<r x - 


- cr 2 ) • <x 3 cr 2 • 


^3^1,! 

2 ' 


= -12, 


(B2) 








- <T 2 ) • A(Tj • 


2 1 


= 0, 


(B3) 






Tr^x 


— <T 2 ) • CT 3 (7j 


■APi 1 

2 1 


= 0, 


(B4) 




Tr (<r 


l — 0*2) ■ 


A (<Ti - (T 2 ) 


2' 1 


= U-B, 

3 


(B5) 


TV 

lr {(Ti 


- cr 2 ; 




— cr 2 J ■ ts <x 3 


■ O ri 1 


- ~y( A x ±»J • C , 


(Bbj 


Tr(<n 




)-A<Jy 


S (<Ti - <T 2 ) 


■ A Pi , 

2 ' x 


= 0, (.7 = 1,2,3), 


(B7) 








Trty (tc) cr fc 


■ A Pi , 

2 


= 0,(i,j,fc = 1,2,3), 


(B8) 








Trt? 2 i 




= 16a; 4 , 


(B9) 








Tr^, 




= 4a: 4 , 


(B10) 


Trt l2 ( 


Cc)t 23 ( 




= Trt 23 (y)t 12 i 




= 4[*V - 3(a • yf] , 


(Bll) 


Trt 23 ( 


30*31 1 




= Trt 3 i(t/)t 23 i 




= -4[a;V - 3(a; • i/) 2 ] . 


(B12) 



2. Proof of the equivalence of redundant states 

First consider states |?/>i) and |t/> 2 ) defined in terms of the two operators in Eq. fl3"0|) by 
= IK-^Xda-pxAxf 2 ^) 

= |[(<ri - <r 2 ) • po-3 • A - (o-a - tr 2 ) • A <x 3 • p] x g 2=1 ) , (B13) 
|-0 2 > = K(<n - <x 2 ) ■ p x A X f 2=1 ), (B14) 

where the spatial wave functions are suppressed for simplicity. Using the trace identities 
shown above, it is easy to show that for j = \ 

16 

<-0i|-0i> = 4(^i> = 2<VaiV2> = 2(^1) = y (p x A) 2 . (B15) 

Thus, if we choose a state such that |0) = \ipi — 2ip 2 ), then ('01V') = an d the states |0i) 
and 2|0 2 ) are equivalent. We will therefore drop the second operator in Eq. (|30"D in writing 
the general baryon wave function. 
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Next, we show that the states given by the action of the operators tu(x) and t23{x) + 
tsi(x) in Eq. ([!!]) on xf 3 12_1 are a l so n °t independent for j — |. Let us consider the following 

state 

|V>> = IM*) +t 23 (x) +t 31 (x)) X % 2=1 ) ■ (B16) 
A calculation of ( , 0|'0) using the trace method yields 

(V#> = ~Tr[t 12 {x) +t 23 {x) +t 31 {x)] 2 Pi x 
z 2 

= ±Tr[t 12 (x)+2t 23 (x)] 2 P hl 

= ]-Ti \t\ 2 (x) + At 2 23 {x) + 2t 12 (cc)t 2 3(a;) + 2t 23 (a;)t 12 (a;)]Pi )1 

L 2 

= ^[16x 2 + 16a; 2 + 2(-8a 2 ) + 2(-8a 2 )] = . (B17) 

Hence, the spin functions ti 2 (x)x^ 2 ~ 1 and [t 23 (x) + t 3 x(x)]xj^ 2 ~ 1 in Eq. (|31| ) are equivalent. 
The same result holds for the operators in Eq. (j3~2"|) . 



APPENDIX C: CALCULATION OF MATRIX ELEMENTS 

Let us consider the spatial matrix element < ^q\V\x/jq ) >, A,B = a, b, c, d, and e. 
Here V = Vb(7y)P(p, A) with Vo is a function of the absolute value of and P(p, A) is a 
polynomial of p and A. The Gaussian wave functions ^{Pa) are defined by Eq. (f45|). Then, 
one has 

B A 3 A B B B B 

(^w\<) = ( ^^ A 7^ A ) 3/2 

x /" dpdAy (r^)P(p,A)exp [-(/3 2 p 2 + ^ 2 A 2 )] , 



ftp 2 + ~ft>? = ft 1 + ^ (P - ™0 a + ^I7« 2 , (C4) 



& 2 ^_ a 2 

2 



a 2 + x6 2 



(CI) 



where 

ft = \(ft 2 + Pf) , ft = \{pf + pf) ■ (C2) 

If Tij = r 12 = p, then the integration is straightforward. 
Now consider the case when 

= a\ + bp , (C3) 
where a, b are constants. If one denotes = s, then A = (s — bp) /a and 
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where x = $\//3p and w = xb/(a 2 + xb 2 ) . A translation p —> p + ws gives the result 
A = (ws + bp) I a with w = 1 — bw. Eq. fld]) now becomes 



?AoAoBoB\ 3 / 2 



7T 2 

x | ^V (|s|)P(p + W s, (*s + 6p)/a)exp [-(/3> 2 + /^s 2 )] 



tf=4 2 (l + 41 , ^ = -^W (C6) 



where 

R' 2 = R 2 (\ \ R' 2 = 

a 2 + xb 2 

At this step, the integration of the right hand side of Eq. (|C5|) becomes straightforward. 

For illustration, let us consider a simple case when ijjfi = ip^ = and V — (p X X) 2 / p' . 
In this case, Vb(ry) = 1/p' and P(p, A) = (p x A) 2 . Since p' = A — p/2, one has a = 1, 
6 = -1/2. Also in this case /3 P = /3 P , /3 A = (3\, x = x = /3f//3f, and w = -2x/(4 + x). Note 
that under the translation p — > p + u;s, the expression (p x A) 2 transforms into (p x s) 2 
that can be replaced by |p 2 s 2 after working out the angular dependent parts. One finds 



2 1 



(C7) 



where 



p? = /?; 2 (i + 1) , = #7 (i + f ) , = aa/v (c? 

We can use similar methods with the momentum-space wave functions 



i \ 3/2 / pI A 



*^* )= toJ ^{-wrw (C9) 

to evaluate integrals which involve the energies Ei = \J pj + mf. Thus, 

U^pi + rn 2 ) = lj (^|) , (l/^/pi+mt) = (^) , (Cll) 



where 



l n (z 2 ) = -^\ dtt 2n+2 VW+^e- t2 , J n (z 2 ) = ^ r dt-^==e- t2 . (C12) 
V 71 " 7o v 71 Jo yt 2 + 2 2 
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The integrals /„ and J n must be evaluated numerically. 

Working out the spin matrix elements and applying the technique of integration shown 
above, one can evaluate the matrix elements (ipi m \H\ipi m ) . For example, the matrix 
element (tJ)q\Ho\i/)q) is given by 

(Vol Wo) = W' P h {m\/(3' p 2 ) + f3Jo (ml/ Pi) (C13) 
a ( 1 2 \ 4a s 1 / <\ 

The matrix elements of Hq in the states in Eq. (p3|) with L p , La 7^ involve the J n with 
n = 0, 1, 2. 
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